Subtraction of test mass angular noise in the LISA Technology Package interferometer 
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We present recent sensitivity measurements of the LISA Technology Package interferometer with 
articulated mirrors as test masses, actuated by piezo-electric transducers. The required longitu- 
dinal displacement resolution of Qpm/vHz above 3mHz has been demonstrated with an angular 
noise that corresponds to the expected in on-orbit operation. The excess noise contribution of 
this test mass jitter onto the sensitive displacement readout was completely subtracted by fitting 
the angular interferometric data streams to the longitudinal displacement measurement. Thus, this 
cross-coupling constitutes no limitation to the required performance of the LISA Technology Package 
interferometry. 



I. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) p^ is 
a joint space mission from the European Space Agency 
(ESA) and the National Aeronautics and Space Adminis- 
tration (NASA), designed as a gravitational wave obser- 
vatory in the frequency range of 0.1 mHz to 0.1 Hz. LISA 
consists of a three spacecraft constellation in a equilat- 
eral triangle formation, flying a total of six free-falling 
test masses that act as end-mirrors of laser interferom- 
eters sensitive to position fluctuations AL better than 
40pm/vIIz over the interspacecraft separation L of 5 
million kilometers. These fluctuations in the separation 
between two test masses are caused by the space-time dis- 
tortion caused by gravitational waves, as well as residual 
acceleration noise. The corresponding strain sensitivity 
AL/L of LISA is of the order of IQ-^i/Vlfe. LISA re- 
quires highly challenging technology that is under devel- 
opment and cannot be tested on Earth. To this end, ESA 
will launch the technology demonstration mission LISA 
Pathfinder (LPF), which consists of a single satellite car- 
rying two payloads: the LISA Technology Package (LTP) 
provided by ESA, and the Disturbance Reduction System 
(DRS) from NASA. LTP is a set of experiments de- 
signed to test core technology essential for LISA, such 
as: 

1. free-fall motion of a test mass with acceleration 



noise lower than 3 x 10 ^*ms ^/vHz at ImHz, 

2. high-precision laser interferometry with a free- 
falling mirror (LTP test mass) with displacement 
sensitivity better than 9 x lO^^^m/v^Hz between 
3 mHz and 30 mHz over a wide dynamic range (sev- 
eral microns), 

3. satellite position correction via micronewton 
thrusters to assure a closed drag-free test mass dis- 
placement control loop. 

4. assess reliability and lifetime of components in 
space, such as optics, and lasers, among others. 
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The main concept of LTP is to shorten one 5 x 10^ m 
LISA interferometer arm to a distance of about 30 cm. 
A laser interferometer is located between the two LTP 
test masses and measures fluctuations in their separation 
with a resolution better than 9 pm/vHz, as well as their 
angular orientation with a sensitivity of lOnrad/vHz. 
The LTP test mass (TM) is a reflecting cube made of 
a platinum-gold (Pt-Au) alloy and resides in a electrode 
housing (EH). The electrodes at the EH internal sides 
and the corresponding faces of the cubic TM form a ca- 
pacitance that can be measured to obtain the TM posi- 
tion. It is also possible to actuate the TM position by 
applying an electric field on the TM through the elec- 
trodes. The Drag-Free Attitude and Control System 
(DFACS) 13 uses the optical metrology output and the 
capacitive sensing as error signals to control a drag-free 
motion of the TM. The micronewton thrusters are the ac- 
tuators on the satellite position to close the DFACS con- 
trol loop. Due to limited gain in the DFACS control loop, 
the test masses have residual angular noise with respect 
to the spacecraft that couples into the longitudinal inter- 
ferometric measurement, thus affecting the performance 
of the optical metrology system. This article presents in- 
vestigations conducted on this effect. As free- falling test 
masses we used articulated mirrors of the interferome- 
ter with 3-axes piezo-electric transducers (PZT) , and the 
engineering model of the LTP optical bench was used as 
optical metrology instrument. 



II. LISA TECHNOLOGY PACKAGE 
INTERFEROMETRY 

The LTP interferometer [4 is divided into two parts: 

• a ultra-stable optical bench, consisting of a set 
of four non-polarizing heterodyne Mach-Zehndcr 
interferometers, whose fused-silica optical compo- 
nents are bonded onto a Zerodur® baseplate [S] , 
what provides high thermal and mechanical stabil- 
ity, and 

• a comparably unstable modulation bench contain- 
ing the laser source and two acousto-optic modu- 



lators that provide two slightly frequency shifted 
laser beams and respective fiber coupling to trans- 
fer the two modulated beams. 

Figure [l] outlines the optical paths of these four interfer- 
ometers, which can be described as follows: 

• The XI 2 interferometer (IFO) measures the fluc- 
tuations in the separation between the two drag- 
free test masses TM1-TM2. Beam 1 and Beam 
2 overlap at the beam combiner BSIO, and the in- 
terference signal is obtained from the redundant 
quadrant photodiodes PD12A and PD12B. 

• The XI IFO monitors the TM 1 position fluctua- 
tions with respect to the optical bench. Both beams 
recombine at the beamsplitter BS8, and the inter- 
ference signal is obtained from the redundant quad- 
rant photodiodes PDIA and PDIB. 

• The reference IFO operates within the ultra-stable 
optical bench only, detecting disturbances common 
to all interferometers that couple into the measure- 
ment in the unstable part (modulation bench and 
fiber optics) , such that they can be subtracted from 
X12 and XI. The recombination beamsplitter for 
this interferometer is BS5 and the readout photode- 
tectors are PDRA and PDRB. 

• The frequency stabilization IFO has an intention- 
ally large optical pathlcngth difference, in order to 
sense the laser frequency noise, and its output sig- 
nal is used to actively stabilize the laser frequency. 
The two beams overlap at BS7 and the input sig- 
nal for the laser frequency control loop is obtained 
from the photodetectors PDFA and PDFB. 

Figure [2] is a photograph of the space-qualified engineer- 
ing model of the LTP optical bench ^J and points out the 
location and mounting of test mirrors that simulate the 
LTP test masses in our experimental setup. The pho- 
tocurrents are processed by a dedicated phasemeter ^ 
that performs a single-bin discrete Fourier transform at 
the heterodyne frequency (difference frequency between 
the two slightly frequency shifted beams) on field pro- 
grammable gate array (FPGA) based digital hardware. 
The phase of each interferometer is computed as the arc 
tangent between two data streams in orthogonal quadra- 
ture. Finally, the longitudinal phase cj) is obtained from 
the difference between the reference phase (phase of the 
reference IFO) and the phase of the measurement inter- 
ferometers (XI, XI 2, and frequency stabilization). The 
conversion from longitudinal phase 4> to TM longitudinal 
displacement ^ is 



* 



4^^ 



(1) 



where A — 1064 nm is the wavelength of the Nd:YAG 
laser used. All photodetectors at the output of all four 



interferometer are quadrant photodiodes (QPD), in or- 
der to obtain alignment information and sense angular 
motion of the test masses by applying a differential wave- 
front sensing (DWS) [51 [S] technique. The horizontal 
wavefront tipping tp can be calculated from the phase 
difference between the left (quadrants A and C) and the 
right side (quadrants B and D) of the QPD. Similarly, the 
vertical wavefront tipping -q is obtained from the phase 
difference between the upper (quadrants A and B) and 
lower side (quadrants C and D) of the QPD. The main 
output of the optical metrology is 

• 'I'l: longitudinal position fluctuations of TM 1 with 
respect to the optical bench, 

• Lfi. horizontal angular motion of TM 1 in the XI 
IFO, 

• rji. vertical angular motion of TM 1 in the XI IFO, 

• ^^12: longitudinal distance fluctuations between 
TMl andTM2, 

• (/9i2: combination of the horizontal angular motion 
of TM 1 and TM2 in the X12 IFO, and 

• 'ni2'- combination of the vertical angular motion of 
TM 1 and TM 2 in the X12 IFO. 

The requirements on the LTP interferometric sensitivity 
have been met with Zerodur® static mirrors [10] shown in 
Figure [2j Diverse noise sources in the system have been 
studied and corrected and the performance has been ex- 
perimentally demonstrated on the engineering model of 
the LTP optical bench [llj . The optical bench is nor- 
mally operated in a vacuum chamber to reduce the effect 
of acoustic noise, and thermal and mechanical fluctua- 
tions in the optical measurement. In order to investigate 
the effect of test mass residual angular noise into the lon- 
gitudinal measurement, the Zerodur® static mirrors were 
substituted by PZT actuated mirrors, which arc shown in 
Figure [3] Figure |4] shows the sensitivity spectra reached 
for the longitudinal TM displacement ^i and ^12 mea- 
sured with Zerodur® static mirrors and with forward bi- 
ased PZT actuated mirrors in static condition. It can 
be seen that all sensitivity curves remain below the re- 
quired 9 pm/-\/Hz in the measurement band (interferom- 
eter budget). Figure [5] shows the horizontal (</?i,i2) and 
vertical (?7i,i2) angular resolution spectra achieved with 
forward biased PZT mirrors in static condition, which is 
better than the required 10 nrad/\/lIz TM jitter in the 
measurement band. 



III. TEST MASS ANGULAR NOISE 
CHARACTERIZATION 

Simulations conducted on the TM dynamics under 
DFACS control led to spectral predictions of the residual 
TM angular noise J12j . From this spectral information, 
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FIG. 1. Optical model of the LTP optical bench engineering model. Note the interference points of the four interferometers. 
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FIG. 2. Space-qualified engineering model of the LTP optical 
bench. Note the expected location of the LTP test masses 
and the mounting of the test mirrors that simulate them. 



FIG. 3. Comparison between PZT actuated mirrors (left) and 
static test mirrors (right). 



we generated a time series that matches this spectral be- 
havior and injected it to the PZT actuated mirrors via a 
digital analog converter (DAC). Figure [6] shows the simu- 
lated angular noise spectra (yi,i2 and ?7i,i2), and the cor- 
responding TM angular noise spectra read out by the in- 
terferometers and computed with a DWS algorithm. Due 
to limited accuracy alignment of the laser beam onto the 
TM center of rotation, the resultant cross-coupling from 
TM angular noise into the TM longitudinal displacement 



readout introduces excessive noise into the measurement 
of TM position fluctuations, thus spoiling the sensitiv- 
ity of the optical readout. Our aim is to characterize 
this cross-coupling and to quantitatively obtain the cor- 
responding coupling factors that translate this angular 
motion of the TM into an apparent longitudinal TM dis- 
placement. Once these coupling factors have been esti- 
mated, the excessive noise can be subtracted from the 
main longitudinal data stream. 
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FIG. 4. Sensitivity spectra of longitudinal phase measure- 
ments performed with static test mirrors and with PZT actu- 
ated mirrors. 
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FIG. 6. Expected injected residual TM angular noise for on- 
orbit operation and interferometric measured angular noise of 
the test masses. 
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5. Sensitivity spectra of angular measurements with 
actuated mirrors. 



A. Fit algorithm 

The longitudinal raw measurements ^i and ^12 de- 
pend on the TM angular motion iy9i,i2 and 771,12 



*i = *i(i,93i,77i) 
*12 = *12 {t,(pi,r]i,(pi2,'r]i2). 



(2) 
(3) 



As it can be seen in Figure |6J the signal-to- noise ratio of 
the TM angular motion is much higher within the LTP 
observation band (3 niHz — 30 niHz) than at higher fre- 
quencies. Under normal laboratory conditions, the effect 
of fast electronic and mechanical noise in the higher fre- 
quency band (approximately above 100 mHz), as well as 
long-term thermal drifts at frequencies below 1 mHz dom- 
inate the time evolution and behavior of the longitudinal 
and angular interferometric signals. Hence, the informa- 
tion of the TM angular noise vanishes in the noise level 
of the measured time series. In order to overcome this 



limitation, we decided to band-pass filter each longitudi- 
nal (^i"^ and ^12) ^"^d angular {ip-^''^,r]-^'^ and </'i2'^i2) 
time series within the LTP observation band. This way, 
we can precisely characterize the cross-coupling process 
of the TM angular motion into the longitudinal interfer- 
ometric readout. The dependency of the filtered longi- 
tudinal measurements ^^p with respect to the filtered 
angular signals (p^^ , r/'^P can be rephrased as 
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Detailed optical simulations have shown a nonlinear cou- 
pling mechanism of parabolic type, but have also shown 
that for the noise levels occuring in our experiment, a 
linear model is sufficient. A general model for this ap- 
proach can be described by the following linear system 
of equations: 

• K = *, (6) 

where is the design matrix for our fitting problem (an- 
gular data f^^, jy'^P), k is a vector containing the coupling 
factors we are looking for, and ^ is a vector correspondig 
to the time series of our target function (longitudinal TM 
data ^'^P). The dimensions of are iV x to, where N 
is the length of the time series and m is the number 
of input time series to be used; in the case of the XI 
IFO TO = 2 (^^'''''71''): and for the X12 IPO to = 4 
((/jj^P, 77]^P, (/'12 1 ^12 )■ '^ is a vector with dimensions m x 1, 
and ^ is a vector with dimensions iV x 1. 
In our specific case, we have the following system of equa- 
tions for the XI IPO: 



*]vx1 = 0]n 
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The system of equations for the X12 IFO can be ex- 
pressed as: 
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(10) 



The fit can be performed by a general hnear least squares 
algorithm. This linear system of equations can be solved 
by applying different algorithms such as the Cholesky de- 
composition, the use of normal equations, or the singular 
value decomposition, among others. The proper selection 
of the solving method usually depends on the topology 
of the design matrix 0. This way, we obtain the set of 
coupling coefficients k^ and k^^ of the TM angular noise 
into the longitudinal TM displacement readout. 
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FIG. 8. Solid curve: sensitivity of the longitudinal phase read- 
out when injecting TM angular noise. Dashed curve: sensi- 
tivity of the corrected longitudinal phase after the angular 
noise subtraction. Dashed curve with crosses: independent 
reference measurement with no injected angular noise. 



B. Test mass angular noise subtraction 



The band-pass filtered data '^^^ and (^'^P,7y^P are uti- 
lized to obtain the coupling coefficients n. For example, 
typical fitted values for k^ are 



[m/rad] 



-4.38 X 10-*^ 
-2.19 X 10-5 



(11) 



Once we have estimated them, it is possible to subtract 
the TM angular noise from the original (unfiltered) mea- 
sured longitudinal TM data as follows: 
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The entire procedure to subtract the TM angular noise 
from the longitudinal interferometric signal is outlined by 
Figure [7] Figure [8] presents the results obtained from this 
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FIG. 7. Flow diagram of the procedure to subtract the TM 
angular noise from the longitudinal phase data stream. 

subtraction. The solid curve is the sensitivity reached by 
the longitudinal phase readout ^ when introducing TM 



angular noise (note that it exceeds the required noise 
budget). The dashed curve for ^I^now is the sensitivity 
achieved after subtracting the fitted angular noise to the 
data of the solid curve. The dashed curve with crosses 
is the sensitivity obtained from an independent measure- 
ment where no angular noise was injected to the test 
masses (PZT actuated mirrors). In general, noise sub- 
traction procedures have to be performed very carefully, 
since there is a non-vanishing probability to corrupt the 
data. In our case, the longitudinal and angular degrees 
of freedom (DOF) are sufficiently orthogonal, such that 
the cross-coupling between them can be very well quan- 
titatively characterized. The linear transformation be- 
tween these two reference systems has been experimen- 
tally measured and can be expressed, for example for 
TMl, as 
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An example of a problematic situation where noise sub- 
traction would be expected to corrupt signal is if ^ cou- 
ples into if, and (p back again into ^P, with factors such 



that. 



to the sensitivity of the interferometric longitudinal TM 
position measurement in LTP. 



9* dip 



(17) 



In our case, however, this product is of the order of 10"^ 
such that no real signal ^ is subtracted. A more detailed 
analysis is currently under investigation. Typical values 
for the TM motion are of the order of 
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As it can be seen in Figure |8j the corrected data reaches 
the same level of the reference measurement where no 
angular noise was applied, which indicates that the com- 
plete cross-coupling effect from the TM angular noise into 
the longitudinal measurement was fitted and extracted 
without corrupting the data. Hence, the residual TM jit- 
ter due to the limited DFACS gain is not a limiting factor 



IV. CONCLUSIONS 

We have presented current sensitivity curves of the 
LTP interferometry, measured at the engineering model 
of the optical bench. We have also demonstrated that 
it is possible to reach picometer resolution in the opti- 
cal readout at a few mHz even with test masses that 
are articulated by forward biased piezo-electric transduc- 
ers. This is an important conclusion for LISA, where the 
test masses will have comparable angular jitter. Further- 
more, we performed experimental investigations on the 
noise contribution of residual test mass angular noise to 
the longitudinal test mass displacement, concluding that 
this cross-coupling process can be fully characterized and 
completely extracted from the longitudinal measurement 
data stream. We obtained coupling factors for the an- 
gular fluctuations, by fitting the measured angular data 
series to the longitudinal data with a linear least squares 
algorithm, using only the expected noise but no addi- 
tional calibration signal. This method can also be used 
in other applications to characterize noise sources of dif- 
ferent kind of systems. 
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